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Abstract 

The deployment of DeepCore array significantly lowers IceCube's energy threshold to about 10 
GeV and enhances the sensitivity of detecting neutrinos from annihilations and decays of light dark 
matter. To match this experimental development, we calculate the track event rate in DeepCore 
array due to neutrino flux produced by annihilations and decays of galactic dark matter. We also 
calculate the background event rate due to atmospheric neutrino flux for evaluating the sensitivity 
of DeepCore array to galactic dark matter signatures. Unlike previous approaches, which set the 
energy threshold for track events at around 50 GeV (this choice avoids the necessity of including 
oscillation effect in the estimation of atmospheric background event rate), we have set the energy 
threshold at 10 GeV to take the full advantage of DeepCore array. We compare our calculated 
sensitivity with those obtained by setting the threshold energy at 50 GeV. We conclude that our 
proposed threshold energy significantly improves the sensitivity of DeepCore array to the dark 
matter signature for < 100 GeV in the annihilation scenario and < 300 GeV in the decay 
scenario. 

PACS numbers: 14.60.Pq, 14.60.St 
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I. INTRODUCTION 



Many astrophysical observations have confirmed the existence of dark matter (DM), which 
contributes to roughly 23% of the energy density of the Universe. Among many proposed 
DM candidates, weakly interacting massive particles (WIMPs) [l, 2| are popular proposals 
since they are theoretically well motivated and also capable of producing the correct relic 
density. WIMPs could annihilate or decay into particles such as electrons, positrons, protons, 
antiprotons, photons, and neutrinos. It is possible to establish the WIMP signature through 
detecting these particles 3-15|. 

Research activities on WIMPs have been boosted recently in efforts of explaining the 



observed anomalous positron excess in the data of PAMELA 
tron excess in the data of FERMI 



and positron plus elec- 



12|. To account for spectral shapes observed by these 



experiments, WIMPs must annihilate or decay mostly into leptons in order to avoid the 
overproduction of antiprotons. This could indicate that DM particles are leptophilic in their 
annihilations or decays [igI [l^]. It has been pointed out that the observation of neutrinos 
can give stringent constraints on the above scenario. Measurements of upward going muons 
by Super-Kamiokande observatory place a limit on the galactic muon neutrino flux, which 
in turn rules out the possibility of WIMP annihilations to r"^r^ as a source of anomalies 



18N20|. Furthermore, one expects that the possibilities of WIMP annihilations into /i^, and 



WIMP decays into fi^ and will all be stringently constrained 2lN23| (see also discussions 



in Ref. 2J]) by the data from IceCube detector augmented with DeepCore array. 



The DeepCore array 25|,|26|] is located in the deep center region of IceCube detector. This 



array consists of 6 densely instrumented strings plus 7 nearest standard IceCube strings. 
The installation of DeepCore array significantly improves the rejection of downward going 
atmospheric muons in IceCube and lowers the threshold energy for detecting muon track 



or cascade events to about 5 GeV. As summarized in Ref. [26|, the low detection threshold 
of DeepCore array is achieved by three improvements over the IceCube detector. First, the 
photo-sensors in DeepCore are more densely instrumented than those of IceCube, as just 
mentioned. Second, the ice surrounding the DeepCore array is on average twice as clear as 
the average ice above 2000 m 27|. Such a property is useful for reconstructing lower- energy 
neutrino events. Finally the DeepCore array uses new type of phototube which has a higher 
quantum efficiency. 
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It is clear that DeepCore array improves the sensitivity as well as enlarges the energy 
window for observing neutrinos from DM annihilations or decays in the galactic halo. Previ- 
ous analyses on the detection of these neutrinos in DeepCore 



23 



28 



have set the threshold 



energy at (40 — 50) GeV for both track and cascade events. For neutrino events with energies 
higher than 50 GeV, the estimation of atmospheric background event rate is straightforward 
since oscillation effects can be neglected. However, to take the full advantage of DeepCore 
array, it is desirable to estimate the track and shower event rates due to atmospheric neu- 
trinos in the energy range 10 GeV < E,^ < 50 GeV. In this energy range, the oscillations 
of atmospheric neutrinos cannot be neglected. In this article, we take into account this 
oscillation effect and calculate the track event rate with a threshold energy E^^ = 10 GeV 
due to atmospheric muon neutrinos from all zenith angles. Given such a background event 
rate, we then evaluate the sensitivities of DeepCore array to the neutrino flux arising from 
DM annihilations and decays in the galactic halo. In the subsequent paper, we shall analyze 
the corresponding sensitivities associated with cascade events. 

This paper will focus on neutrino signature induced by low-mass DM. Hence our interested 
DM mass range is far below TeV level implied by PAMELA and FERMI data. Therefore 
we shall consider neutrino flux induced by DM annihilations/decays into both leptons and 
hadrons. Specifically, we consider the channels XX ~^ bb, t~^t~ , and fi~^fi~ for annihilations 
and the channels x ~^ bb, t^t~ and for decays. Since we are only interested in low- 

mass dark matter, we have neglected neutrino fluxes generated through DM annihilations 
or decays into tt, W^W~ and ZZ final states. We also neglect neutrino fluxes arising from 
light meson decays, as the annihilation cross section for xX ~^ is likely to be suppressed 
by rrig j^. We shall compare the constraints on DM annihilation cross section and DM 
decay time for different values of threshold energy Ejj^^. For such a comparison, we employ 
the modes xX ~^ fi~^fi~ and x fi~^fi~ for illustrations. 

This paper is organized as follows. In Sec. II, we outline the calculation of muon neutrino 
flux from WIMP annihilations and decays in the galactic halo. In Sec. HI, we calculate the 
atmospheric muon neutrino flux from all zenith angles with E^, > 10 GeV. The oscillations 
between z/^ and u^- are taken into account. In Sec. IV, we evaluate the sensitivity of DeepCore 
array to neutrino flux arising from WIMP annihilations or decays in the galactic halo. We 
compare our results with those obtained by setting E^^ = 50 GeV. We summarize in Sec. 
V. 
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II. NEUTRINO FLUX FROM ANNIHILATIONS AND DECAYS OF DARK MAT- 
TER IN THE GALACTIC HALO 



The differential neutrino ffux from the galactic dark matter halo for neutrino flavor i can 



be written as 



20| 



d$,, An (av) 



dE„ 



4tt 2m2 



for the case of annihilating DM, and 
d$,.- An 1 



dE 



dE 



(1) 



RqPq X Ji(Af]) 



(2) 



for the case of decaying DM, where Rq = 8.5 kpc is the distance from the galactic center 
(GC) to the solar system, pQ = 0.3 GeV/cm'^ is the DM density in the solar neighborhood, 
is the DM mass, is the DM decay time and dN^JdE is the neutrino spectrum 



per annihilation or decay for a given annihilation or decay channel F with a corresponding 
branching fraction Bp- The neutrino spectra dNj^JdE for different channels are summarized 



in Ref 



28 



|29|. The quantity {(tv) is the thermally averaged annihilation cross section, which 
can be written as 

{av)=B{av)o, (3) 



with a boost factor B 30, ImI]- We set {(yv)^ = 3 x 10~^^ cm^s~^, which is the typical 
annihilation cross section for the present dark matter abundance under the standard thermal 
relic scenario We treat the boost factor i? as a phenomenological parameter. The 
dimensionless quantity J„(Af2) is the DM distribution integrated over the line-of-sight (l.o.s) 



and averaged over a solid angle Af2 = 2tt{1 — cosipj 

1 



I.e., 



An 



An 



l.o.s ^& 



Pq 



(4) 



where p is the DM density at a specific location described by the coordinate (/, ip) with / the 
distance from the Earth to DM and ip the direction of DM viewed from the Earth with ip = 
corresponding to the direction of GC. The distance r = a/ Rq + P — 2RqIcos'iP is the distance 
from GC to DM. The upper limit of the integration, /^ax = RqCos^Ij + \ Rj. — -R|sin^?/;, 
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depends on Rg, the adopted size of the galactic halo. In this analysis, we take Rg = 20 kpc 



and use the Navarro-Frenk- White (NFW) DM density profile 

2 



p(r) 



Rg 



Rg 



Rg + r 



(5) 



with ps = 0.26 GeVcm"^ such that = 0.3 GeVcm"^. 

Neutrinos are significantly mixed through oscillations when they travel a vast distance 
across the galaxy. We determine neutrino flavor oscillation probabilities in the tribimaximal 
limit 33| of neutrino mixing angles, i.e., sin^ ^23 = 1/2, sin^0i2 = 1/3 and sin^ ^13 = 0. The 



neutrino fluxes on Earth are related to those at the source through 



34- 



$^ = ^$0 ^ 2 2 



and 

2 7 7 

9 18 U 

where is the neutrino flux of flavor i at the astrophysical source. It is understood that 



= '^u. = + ^< + (7) 



the recent T2K [37[ and Double Chooz [38[ experiments have indicated a non-zero value 
for 613. Taking the T2K best-fit value sin^ 26'i3 = 0.11 at the CP phase 6 = for the normal 
mass hierarchy, we have 



.0 , n 0^7-1)0^ 0.37$°^ 
^0 I n /io,T>0 



$ = 0.53$" +0.26$" +0.21$^ , 

•^e Ufi //// t^T ' 



= 0.26$" + 0.37$" + 0.37$^ , 

"ii i^a ' 



$^ = 0.21$" + 0.37$" + 0.42$" . (8) 

To proceed our discussions, let us first take the neutrinos at the source to be those generated 
by B meson decays following the xx ~^ bb annihilation. In this special case, $"^ = $"^ = $"^ 
at the source, and consequently the relation $1,^ = $1,^ = $^^ always holds due to the 
probability conservation, irrespective of the form of oscillation probability matrix. Let us 
now consider neutrinos produced at the source by muon decays following the xX ~^ 
annihilation. In this case, one has $" = $" and $" = 0. Taking $" = $" = $", 
one obtains = 0.78$" and $,^^ = $,^^ = 0.61$" for tribimaximal values of neutrino 
mixing parameters. On the other hand, with the T2K best-fit 6*13 value, one arrives at 
$1.^ = 0.79$", $^ = 0.63$" and = 0.58$" for the normal mass hierarchy. Clearly $y^ 



is almost unaffected while (^y^/^y^ — 1 = 9%. For the inverted mass hierarchy, one obtains 
the same z/g ffux while ^v^^/^vr ~ 1 = 12%. Hence T2K result implies an 0(10%) difference 
between the arrival and Vr ffuxes for neutrinos produced by XX ~^ fJ^^fJ^" annihilations. 



III. ATMOSPHERIC NEUTRINO FLUXES 



Knowing atmospheric neutrino background is important for evaluating the sensitivity of 
DeepCore array to neutrino flux from DM annihilations or decays. We begin by computing 
the flux of intrinsic atmospheric muon neutrinos arising from pion and kaon decays, following 
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the approaches in Refs 

dEdX 



40|. The z/^ flux arising from vr decays can be written as 



r^N Q{E^ - ^) 
E Je d^E^[l - 7^) 

— PiriEn, X, X')—Fi^^{E^, En) 



N 



E^ 



X' 



X exp ( -— ) (j)N{EN), 



(9) 



where E is the neutrino energy, X is the slant depth in units of g/cm^, ^ is the zenith 
angle in the direction of incident cosmic-ray nucleons, = m'^/m'^, d^^ is the pion decay 
length in units of g/cm^, Aa? is the nucleon interaction length and Aat is the corresponding 
nucleon attenuation length. The function P^^lEj^, X, X') is the probability that a charged 
pion produced at the slant depth X' survives to the depth X (> X'), which is given by |43| 

where A^r = 160 g/cm^ is the pion attenuation length, r^r is the pion lifetime at its rest frame. 



while p{T) is the atmosphere mass density at the slant depth T. Finally, Fnt^ 
the normalized inclusive cross section for N + air ^ n'^ + Y, which is given by 

Et^ d(T(i?^, Ejy^ 



Er,, Em) is 
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ENn{E.„, E 



N 



cM-xr+ + cM 



X 



(11) 



dE^ 

with X = E^/En, c+ = 0.92, c_ = 0.81, p+ = 4.1, and p_ = 4.8. 

The primary cosmic-ray spectrum (PnIEn) in Eq. iQ includes contributions from cosmic 
ray protons and those from heavier nuclei. We have (Pn^En) = J^A^'^^i^'^) with A the 
atomic number of each nucleus. The spectrum of each cosmic-ray component is parametrized 
by 



^AiE^) = Kx{En + 6exp[-c/E^])' 



(12) 



TABLE I: Parameters for all five components in the fit of Eq. (|12p . 



parameter / component 



a 



K 



Hydrogen (A=l) (< lO^GeV) 
Hydrogen (A=l) (> lO^GeV) 
He (A=4) 
CNO (A=14) 
Mg-Si (A=25) 
Iron (A=56) 



2.74 
2.71 
2.64 
2.60 
2.79 
2.68 



14900 
14900 
600 
33.2 
34.2 
4.45 



2.15 
2.15 
1.25 
0.97 
2.14 
3.07 



0.21 
0.21 
0.14 
0.01 
0.01 
0.41 



in units of m ^sr ^GeV ^. T 
They are tabulated in Table I 



le fitting parameters a, K, b, c depend on the type of nucleus. 



42| . The kaon contribution to the atmospheric z/^ flux has 



the same form as Eq. (Q with an inclusion of the branching ratio B{K — > fiu) = 0.635 and 
appropriate replacements in kinematic factors and the normalized inclusive cross section. 
In particular, Fmk{Ek, E^) can be parametrized as Eq. ( fTTl) with c+ = 0.037, c_ = 0.045, 
p+ = 0.87, and p_ = 3.5. 

Since our interested energy range is as low as 10 GeV, the three-body muon decay con- 
tribution to the atmospheric flux is not negligible, particularly in the near horizontal 
direction. To obtain this part of contribution, we first compute the atmospheric muon flux 
from pion and kaon decays. The muon flux induced by pion decays is given by jsol, l4o| 



dE 



X 



poo rt^N 

/ dEM / dE^ / dX"P^{E,X,X") 
Je' Je' Jo 



(i^i?^(l 



'ML 

'IZL. 



E^ 



X / — — P^(i?^,X ,X')—Fn^{E^,En) 
An 







X exp 



E„ 



X^ 

An 



>n{En 



(13) 



where E' and E are muon energies at slant depths X" and X respectively, while Pi_i{E, X, X") 
is the muon survival probability. The muon flux induced by kaon decays can be calculated 
in a similar way. Since fi~{fi~^) produced by 7r^(7r"'") decays are polarized, we classify muon 



flux into four different components sucli as dN^+ /dE, dN^_ /dE, dN^+ /dE, and dA^'^_ /dE. 

f^R t^R /^L 

We also calculate additional four components of the muon flux arising from the kaon decays. 



Hence the z/^ flux arising from muon decays can be written as 40|, |43| 



dEdX ^^Je ^ d,{E^,X)E, 

(14) 



div^±(E^,e,x) 



d^^ 



where d^{E^,X) is the muon decay length in units of g/cm at the slant depth X and 
E^±_^j^ (E/E^) is the normalized decay spectrum of fif — ?■ u^. Summing the two-body and 



three-body decay contributions, we obtain the total intrinsic atmospheric muon neutrino 
flux. In Fig. 1, we show the comparison of angle- averaged atmospheric muon neutrino flux 



obtained by our calculation and that obtained by Honda et al. 



4^. At = 10 GeV, two 



calculations only differ by 3%. At E,y = 100 GeV, the difference is 10%. We also show in the 



same figure the atmospheric muon neutrino flux measured by AMANDA-II detector 45|]. It 
is seen that both calculations agree well with AMANDA results. 

To completely determine the atmospheric muon neutrino flux, one also needs to calculate 
the intrinsic atmospheric tau neutrino flux, although this part of contribution is rather small. 
The intrinsic atmospheric z/^ flux arises from decays. This flux can be obtained by solving 



cascade equations [40|, |46| . We obtain 



d^N,^iE,X) _ Znd^Zd^,^ exp(-X/Ajv)0;v(Ejv) 
dEdX 1-Znn ' An ' 

where Znn = 1 — \n/ An and ZnDs is a special case of the generic expression 

Z.,{E,) = J^^ d^.;^^ 51^ . (16) 

with driiA^jYiEi, Ej) = daiA^jY{.Ei, Ej) / ai^i^Ei) and Aj the interaction length of particle i 
in units of g/cm^. The decay moment Zo^v^ is given by 



oo 



Zd.uAEv^) = / dEo, '^f" \ J?"} r^'^j/" \ Fds^v^ {Ev^ I Ed, ) , (17) 

Je., <PDAEur) dDs[EDj 

where is the decay length of Ds and Fo^^^^^Ey^ / E^^) is the normalized decay distribu- 
tion. In this work, we employ the next-to-leading order (NLO) perturbative QCD with 
CTEQ6 parton distribution functions to calculate the differential cross section of NA — cc 
and determine ZnDs- 
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Fig. 1: The comparison of angle-averaged atmospheric muon neutrino (i/^ + z^^) flux obtained 



by our calculation and that obtained by Honda et al. 
AMANDA-II measurements [45|] is also shown. 



441 ]. Angle-averaged + flux from 



Finally, the atmospheric Vn flux taking into account the neutrino oscillation effect is given 



by 



dE 



dX 

d^K 



dEdX 



■P. 



dEdX 



(18) 



where P,^^,^ (E, L{X, 0) = Pu^-^u, {E, L{X, 0) = sin^ 2^23 sm\l.27Aml^L/ E) is the 
Vt- oscillation probability and L(X, ^) is the linear distance from the neutrino production 
point to the position of IceCube DeepCore array. The unit of Amg^^ is eV^ while L and E are 
in units of km and GeV respectively. The best-fit values for oscillation parameters obtained 



from a recent analysis 48|| are Amgj^ = 2.47 ■ 10^3 and sin 2^23 = 1 respectively. 
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IV. RESULTS 



In IceCube DeepCore, the event rate for contained muons is given by 

= / dE^ dE.^NAPiceVtr 

where pice = 0.9gcm~^ is the density of ice, = 6.022 x 10^'^ is Avogadro's number, 



Vtr ~ 0.04 km'^ is the effective volume of IceCube DeepCore array for muon track events 25|. 
d^u^/dE^^ is the muon neutrino flux arrived at IceCube, -Emax is taken as for annihilation 
and ?Ti^/2 for decay, and E^J^ is the threshold energy for muon track events. In this work. 



we use differential cross sections da^^{Ey^, E^) / dE^j_ given by Ref. |49[ with CTEQ6 parton 



distribution functions. We also set E^}^ = 10 GeV. 

As stated before, we consider neutrino fluxes generated by the annihilation channels 
XX ~^ bb, T~^T", and fi~^fi~, and the decay channels x ~^ bb, t^t^ and Given 
the atmospheric neutrino background, we present in Fig. 2 the required DM annihilation 
cross section as a function of for threshold energy E^^ = 10 GeV and a cone half- angle 
■^max = 1° such that the neutrino signature from DM annihilations can be detected at the 
2(7 significance in five years. Non-detection of such a signature would then exclude the 
parameter region above the curve at the 2(T level. We have presented results corresponding 
to different annihilation channels. One can see that the required annihilation cross section 
for 2cr detection significance is smallest for XX ~^ fJ'^fi^ channel and largest for the channel 
XX ~^ bb. We also present the 3a constraint on xX ~^ t^t~ annihilation cross section 
obtained from Super-Kamiokande data of upward going muons |l^, which has been used 



to rule out W 



MP annihilations into r+r as a possible source of previously mentioned e 



± 



anomalies [18|-|20|. Such a constraint can be compared with the expected 2a constraint on 
the same annihilation channel from DeepCore detector. 

Constraints on DM annihilation cross section were also obtained from gamma ray ob- 
servations and cosmology. H.E.S.S. telescope performed a search for the very-high-energy 
(> 100 GeVj_7-ray signal from DM annihilations over a circular region of radius 1° centered 



at the GC 50|]. With DM particles assumed to annihilate into qq pairs, the limit on DM 



annihi! 
Ref. 



ation cross section as a function of rriy for NFW DM density profile is derived in 



50|. We present this constraint in Fig. 2 as well. For > 300 GeV, the parame- 
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Fig. 2: The dashed hne, thin sohd Une, and dot-dashed hnes are the expected constraints on 
DM annihilation cross section by DeepCore detector for XX ~^ bb, XX ~^ fJ-^f^'^ and XX ~^ t'^t~ 

le H.E.S.S constraint on the annihilation cross 



channels, respectively. The thick solid line is tl 
section of DM into the light quark pair xX ~^ QQ 

XX ~^ e~^e~ annihilation cross section from the analysis of cosmic microwave background (CMB) 



501 ]. The dot-dot-dashed line is the constraint on 



data 



5l[ |. The dotted line is the 3a constraint on the annihilation cross section of xX ~^ t^"^ 



from Super-Kamiokande data 



ter space with B > 100 (i.e., (cry) > 3 x lO^^'^m'^s"^) in Fig. 2 could be excluded by the 
H.E.S.S. data. However, the H.E.S.S. constraint on xX ~^ QQ becomes much weaker for 
< 300 GeV. We point out that this constraint is obtained with NEW profile normahzed 
at Pq = 0.39 GeVcm""^. The H.E.S.S. constraint would be slightly less stringent if our 
adopted normalization Pq = 0.3 GeVcm"^ is used. 

Cosmological Constraints on DM annihilation cross section can be obtained from the data 
of big-bang nucleosynthesis (BBN) and cosmic microwave background (CMB). In such an 
analysis, the DM annihilation cross section is assumed to be velocity dependent such that 
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Fig. 3: The required DM annihilation cross section (xx ~^ A*"*^^") &s a function of such that 
the neutrino signature from DM annihilations can be detected at the 2a significance in five years. 
Results corresponding to different V'maf presented. For comparison, we also show the result 
with Ej}" = 50 GeV and ^Amax = 10° 



5l| 



{-V) = -^Tl^, (20) 
where Vq is the DM velocity at the freeze-out temperature, while the values for e and n depend 



on specific models. For Sommerfeld enhancement [30| of the DM annihilation cross section 
induced by light-scalar exchange, one has n = 1 and e ~ m^/m^ with the light-scalar 
mass. The CMB anisotropy can be affected by the energy injection in the recombination 
epoch due to DM annihilation process such as xX ~^ e~^e~ and xX ~^ W~^W~ . In Fig. 2, we 
show the upper bound on (av) for xX ~^ e~^e~ channel for n = 1 and Tkd = 1 MeV with Tkd 
the kinetic decoupling temperature. This upper bound is inferred from the upper bound on 
crv)o such that the resulting CMB power spectrum remains consistent with observations 



511 . The above upper bound on (erf )o is shown to be sensitive to the parameter e while the 



corresponding bound on (av) is insensitive to it. It will be interesting to convert the above 
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bound on e^e~)v) into the one on {(t{xx ~^ fi'^fi~)v). However such a conversion 

is highly model dependent which is beyond the scope of the current work. 

Having compared the expected sensitivities of DeepCore detector with other experimental 
constraints on various DM annihilation channels, we discuss how the DeepCore constraint on 
DM annihilation cross section varies with the chosen cone half-angle and threshold energy. 
We use the channel XX ~^ fJ^'^fi^ to illustrate these effects. Fig. 3 shows the required DM 
annihilation cross section {(t{xX ~^ At"*"/!")!;) for a 2cr detection in five years for different 
cone half-angle ipma,x- One can see that the constraint on the DM annihilation cross section 
gets stronger as ipniax increases from 1° to 2°. However, the constraint turns weaker as ipuiax 
increases further. This is due to the factor J2(^^)^^ which depends on the square of DM 
density (see Eq. (jl])). The constraint curve rises with an increasing V'max for ipraa.x > 2°, 
since the signal increases slower than the background does for such a V'max range. In this 
figure, we also show the result for a higher threshold energy E]^ = 50 GeV with a cone 



half-angle ipraa,x = 10° for comparison. This result is taken from Ref. [28[ where ipmax = 10° 
is identified as the most optimal cone half-angle for constraining DM annihilation cross 
section at that threshold energy. We note that, for large m^, lowering from 50 GeV 
to 10 GeV results in more enhancement on the event rate of atmospheric background than 
that of DM annihilation. Hence, the constraint on DM annihilation cross section is weaker 
by choosing E^J^ = 10 GeV. On the other hand, for small m^, lowering i?*^ enhances more 
on the event rate of DM annihilations than that of atmospheric background. For < 100 
GeV, one can see that the constraint on DM annihilation cross section with E^}^ = 10 GeV 
is always stronger than that with E^^ = 50 GeV. We note that DeepCore constraints on 
other annihilation channels have similar cone half-angle and threshold energy dependencies. 

Besides studying DeepCore constraints on DM annihilation channels, we also present 
constraints on DM decay time for x ~^ bb, t^t~ and fi~^fi~ channels. Fig. 4 shows the 
required DM decay time for a 2a detection of neutrino signature in five years for each 
channel. We have taken E^^ = 10 GeV and tpmax = 90°. Non-detection of such a signature 
would then exclude the parameter region below the curve at the 2cr level. For comparison, we 
also show 3(7 limit on x ^ t^t^ from Super-Kamiokande data of upward going muons 18|. 
One can see that the channel x ~^ fi'^fJ'~ requires the smallest decay width to reach the 2a 
detection significance in five years of DeepCore data taking. 

Finally we present how the DeepCore constraint on DM decay time varies with the 
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m [GeV] 



10^ 



Fig. 4: The dot-dashed hne, sohd hne and doted Une are the required DM decay time for a 2a 



detection of neutrino signature in five years for x ~^ bb, fj,^ i^i anc 
The dashed line is the Super-Kamiokande constraint on x — ?> T"^r" 



T~^T channels, respectively. 



chosen cone half-angle and threshold energy. We use the channel x ~^ fJ^^fJ^^ to illustrate 
these effects. Fig. 5 shows the required DM decay time {x — fJ'~^f^~) as a function of 
for different cone half-angle V'max such that the neutrino signature from DM decays can be 
detected at the 2a significance in five years. For DM decays, the curve rises with increasing 
V'max since the event rate of DM signal increases faster than that of atmospheric background 
as iprnax lucreases. For comparison, we show the required DM decay time for a 2a detection 
in five years with E^^ = 50 GeV and V^max = 50°. It has been pointed out in Ref. 28|] that 
V'max = 50° gives the most stringent constraint on DM decay time for E^^ = 50 GeV. One 



can see that the constraint on DM decay time is strengthen by lowering E^^ from 50 GeV 



to 10 GeV for < 300 GeV. 
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Fig. 5: The required DM decay time (x — > ^^~^^^~) as a function of such that the neutrino 
signature from DM decays can be detected at the 2a significance in five years. Results corresponding 
to different '^Wx are presented. For comparison, we also show the result with E}j^ = 50 GeV and 



Ipma.^ = 50' 



V. SUMMARY 

We have calculated the track event rate in IceCube DeepCore array resulting from muon 
neutrino flux produced by annihilations and decays of dark matter in the galactic halo. 
In this calculation, we have employed NEW profile for dark matter mass distribution and 
consider the channels XX ~^ bb, t'^t~, and fi^fi~ for annihilations and the channels x ~^ 
bb, T^T~ and fi^fi" for decays. We also calculated the track event rate due to atmospheric 
background. We compare the signal event rate with that of the background for E^^ > 10 
GeV. 

We have presented sensitivities of IceCube DeepCore array to neutrino flux arising from 
dark matter annihilations and decays. For a given dark matter mass, we evaluated the 
dark matter annihilation cross section and dark matter decay time such that a 2a detection 
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significance for tlie above signatures can be achieved by DeepCore array for a five-year data 
taking. Tlie DeepCore sensitivities on dark matter anniliilation cross section were compared 
witli the constraint obtained from H.E.S.S. gamma ray observations and the constraint 
derived from the data of CMB power spectrum. Using XX ~^ and x ~^ 

examples, we also presented how DeepCore constraints on dark matter annihilation cross 
section and dark matter decay time vary with the chosen cone half-angle and threshold 
energy. We like to point out that our calculated sensitivities based upon E^^ = 10 GeV are 
significantly more stringent than those obtained by taking E^J^^ = 50 GeV for m^^ < 100 GeV 
in the annihilation channel and < 300 GeV in the decay channel. 
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